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LOCALITY OF THE THOMAS-FERMI-VON WEIZSACKER 

EQUATIONS 

F. Q. NAZAR AND C. ORTNER 


Abstract. We establish a pointwise stability estimate for the Thomas-Fermi-von Weiz- 
sacker (TFW) model, which demonstrates that a local perturbation of a nuclear arrange¬ 
ment results also in a local response in the electron density and electrostatic potential. 
The proof adapts the arguments for existence and uniqueness of solutions to the TFW 
equations in the thermodynamic limit by Catto et al. (1998). 

To demonstrate the utility of this combined locality and stability result we derive sev¬ 
eral consequences, including an exponential convergence rate for the thermodynamic 
limit, partition of total energy into exponentially localised site energies (and conse¬ 
quently, exponential locality of forces), and generalised and strengthened results on the 
charge neutrality of local defects. 


1. Introduction 

Locality properties of electronic structure models are a key premise in certain state of 
the art numerical algorithms. A well-established example is “near-sightedness”, a locality 
property of the density matrix which gives rise to linear scaling algorithms for Kohn-Sham 
type models [25, 38, 24, 6]. A stronger notion is the locality of the mechanical response, 
which is a fundamental premise underpinning the construction of interatomic potentials 
and of multi-scale algorithms such as hybrid QM/MM schemes [17] (here, it is termed 
“strong locality”). This latter category of locality is less well studied, the only result in 
this direction being the locality of non-selfconsistent tight binding models [14]. 

The aim of the present work is to establish the locality properties satished by the 
Thomas-Fermi-von Weizsacker (TFW) model. Our main technical result to achieve this 
is the following pointwise stability estimate for the TFW equations, which establishes 
the locality of the electron response to changes in the nuclear conhguration. Compared 
with [14] it is noteworthy that our result takes Coulomb interaction fully into account. 
Rigorous statements, under different conditions, are given in Theorems 3.4 and 3.5. 

Theorem. Fori = 1,2 let rrii G represent nuclear charge distributions satisfying 

1 f 

ruj > 0 and lim -- inf / mAz) dz = -|-oo. 

H-i>oo R xGRS 

Let the corresponding ground state electron densities and electrostatic potentials, denoted 
by Ui, —)■ M, satisfy the TFW eguations, 

-Aui + - AiUi = 0, 

-A0i = 47r(mj - uf). 
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Then there exists C, 7 > 0 such that for all y 

\{ui - U2)iy) \ + 1(01 - 02)( 2 /)I < c J{mi - 


( 1 . 1 ) 


In the remainder of the article we explore some of the consequences of this locality result: 
In Proposition 4.1 we obtain new estimates on hnite-domain approximations which yield 
exponential decay of surface energies as well as an exponential convergence rate for the 
thermodynamic limit. In Corollary 4.2 we show that (1.1) gives rise to rigorous results 
that match, and substantially generalise, the Thomas-Fermi theory of impurity screening 
in metals [2, 28]. In Theorem 4.3 we strengthen existing results on the neutrality of 
the TFW model [10]. In all these results, general (condensed) nuclear arrangements are 
treated. 

A striking application of (1.1) is that it allows us to decompose energy into local con¬ 
tributions from which we obtain local site energy potentials: Given a countable collection 
of nuclei Y = (yj)jgi^ C we construct an energy density S{Y-,x) which allows us to 
dehne the TFW energy f^£(Y; x) dx of an arbitrary volume C in a meaningful way. 
This then motivates us to dehne site energies 


Ej{Y) := f ipj{x)S{Y]x) dx, 

Jr3 


where {(pj)j^-^ is a smooth partition of unity of which can be constructed in such a 
way that Ej are permutation and isometry invariant and most crucially, Ej are local in 
the sense that 


aEj(Y) 




( 1 . 2 ) 


for some C, 7 > 0. The rigorous statement of this result is given in Theorem 4.4. An 
analogous result has recently been proven for a tight binding model in [14]. 

This result not only gives a strong justihcation for the construction of classical short- 
ranged interatomic potentials in metals, but in fact it allows us to treat the TFW me¬ 
chanical response as if it emanated from such a classical potential. For example, (i) the 
analysis of the Cauchy-Born continuum limit [37] applies directly to the TFW model; 
and (ii) we can generalise in [13] the analysis of variational problems for the mechanical 
response to defects in an inhnite crystal [19]. 

The remainder of this article is organised as follows: In Section 2 we recall the dehnition 
of the TFW model and summarise the relevant existing results. In Section 3 we state the 
main technical results, including the rigorous statement of the stability result (1.1). In 
Section 4 we present applications. Concluding remarks are made in Section 5, followed 
by the detailed proofs of the results in Section 6 . 


Remark 1. We conclude the introduction with a remark about the analytical context 
of this work. The TFW equations for the electron density and potential is a coupled 
Schrodinger-Poisson system. Other systems of this class can be found in semiconductor 
physics [34, 43, 45]. We also note that Thomas-Combes estimates give conditions under 
which eigenfunctions of a Schrodinger operator decay exponentially [1, 16]. While the 
results obtained for these systems are similar, the corresponding equations have different 
structure, hence the analytical techniques used to study them differ considerably. 

The closest existing result to (1.1) we have found is [ 8 , Theorem 4.6], which shows the 
exponential decay of the electron density away from the boundary of a crystal. Both 
(1.1) and [ 8 , Theorem 4.6] utilise the uniqueness of the TFW equations to prove stability 
estimates. In Section 4.1, we use Proposition 4.1 to generalise [ 8 , Theorem 4.6]. 
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2. The TEW Model 
For p G [1, cxo] we define the fnnction spaces 

;= { / : m3 ^ M I V a C compact, / G A(A) } and 

:= { / e I sup \\f\\LP{B^{x)) < oo }. 


For fc G N, A]^^(M3), are dehned analogonsly. For a mnlti-index a = (oi, 02, 03), 

we define the partial derivative d°‘ = Thronghont this paper, a, (3 denote 

three-dimensional mnlti-indices. 

The Conlomb interaction, for f,g G L®/3(^m 3), is given by 

D{f,g)= f f 

Jrs Jrs \x — y\ JR 3 V I ly 

This is finite dne to the Hardy-Littlewood-Sobolev estimate [3] 

\D{f,g)\ < C||/||2,6/5(R3)||5'||x,6/5(K3). (2.2) 

Let m G L®/3(^m 3) , m > 0 denote the charge density of a hnite nnclear clnster, then the 
corresponding TFW energy fnnctional is dehned, for v G by 

E™(n,m) = C'w/ \Vv\^ + Ctp [ + lD{m - v^,m - v^). (2.3) 

Jrs Jr3 2 

The fnnction v corresponds to the positive sqnare root of the electron density. The hrst 
two terms of m) model the kinetic energy of the electrons while the third term 

models the Conlomb energy. We remark that this dehnition of the Conlomb energy is 
only valid for smeared nnclei. We can rescale the energy to ensnre Cw = Ctf = 1- 

To constrnct the electronic gronnd state for an inhnite arrangement of nnclei (e.g., crys¬ 
tals), we restrict admissible nnclear charge densities to m G LWf(M3),m > 0, satisfying 

(HI) snp / m(z) dz < 00 , 

xGR3 JBi (x) 

1 r 

(H2) lim inf — / m{z) dz = 00 . 

R-^ooxeR^ R JBii{x) 

The property (HI) gnarantees that no clnstering of inhnitely many nnclei occnrs at any 
point in space whereas (H2) ensnres that there are no large regions that are devoid of 
nnclei. 

Let und dehne the trnncated nnclear distribntion = m ■ Xbb.^(o)i then the 

minimisation problem 


rTFW/ 


(m^J = inf <J A™(n,mij„) 


V G iL3(M3),n > 0, 


v^ = 




possesses a nniqne minimiser The charge constraint ensnres that the system is 

nentral. Farther, solves the corresponding Enler-Lagrange eqnation, which can be 
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expressed as the coupled system 


-Amk„ + - f)RnUR„ = 0, 

(2.4a) 

= 47r(mR^ -M^J. 

(2.4b) 

By the proof of [12, Corollary 2.7, Theorem 6.10], it follows that 


l|wRj|Hi„if(R3) + ||0i?„||L=„.j(R3) < C, 

(2.5) 


where C is independent of Rn- Consequently, (2.5) implies that along a subsequence 
{uR^,(f)R^) converge to {u,(l)). Passing to the limit in (2.4) yields the following result. 


Theorem 2.1 ([12, Theorem 6.10]). Let m G L^jjj£(M^),m > 0 satisfy (H1)-(H2), then 
there exists a unique solution (m, 0 ) G L°°(M^) x L^^££(M^), up to the sign of u, of 

5 

—Am + — 4>u = 0, (2.6a) 

o 

—A(/) = 47r(m — M^), (2.6b) 

in the distributional sense. In addition, infw > 0. 

Definition 1. For any nuclear configuration m satisfying (H1)~(H2), we refer to 
solving (2.6) as the ground state corresponding to m. □ 

Remark 2. A concise proof of uniqueness of the TFW equations is given in [7] under 
the assumption that m is smooth and hence u,(j) ^ 1T^’°°(M^), which simplihes the earlier 
proof given in [12]. □ 

In the next section, we discuss results that can be obtained by generalising the proof 
of Theorem 2.1. 


3. Main Results 


3.1. Uniform regularity estimates. In the proof of our main results in the next section 
we employ regularity estimates that rehne those of [12], and may be of independent 
interest. 

Other than Proposition 3.1, our estimates rely on uniform variants of (H1)-(H2) and 
it turns out that (H2) may then be simplihed without loss of generality; see Lemma 6.1 
for more details. Given M,uo,ui > 0, let u = {uo,ui) and dehne the class of nuclear 
conhgurations 


A4i.(M,u;)= mGLLf(®') 




< M , 


VR>0 inf 

xGR® 


Br{x) 


m{z) dz > UqR^ — wi >. (3.1) 


As each nuclear distribution m G Aii, 2 (M,uj) satishes (H1)-(H2), Theorem 2.1 guaran¬ 
tees the existence of corresponding ground states (m, 0). We adapt the proof of existence of 
Theorem 2.1 to show that the uniformity in upper and lower bounds on m G M.l 2 [M,u:) 
yields regularity estimates and lower bounds on these ground states which are also uni¬ 
form. 


Proposition 3.1. For any nuclear distribution m : ^ K> 0 j satisfying 
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there exists {u, (j)) solving ( 2 . 6 ) and satisfying u> 0 and 

^ ^*(1 + (3.2) 

ll'^llH2^if(K3) + (3.3) 

Proposition 3.2. There exists cm,u} > 0 such that for all m G Ml 2 [M,u}) the corre¬ 
sponding ground state {u, (f) is unigue and the electron density u satisfies 

inf u{x) > Cm, to > 0. (3.4) 


Assuming higher regularity of the nuclear distributions implies higher regularity of the 
ground state. We therefore dehne, for A; G Mq, 


J^ljk (iff, cu) 


m e 



< M, 


V i? > 0 inf / miz) dz > ujqB? — cui i. 

^ 6 ®® Jbr{x) j 

Arguing by induction and applying the uniform lower bound (3.4) yields the following 
result. 


Corollary 3.3. Suppose k E Nq and m G , then the corresponding solution 

{u, 0 ) to ( 2 . 6 ) satisfies 

<C{k,M,uj). (3.5) 


3.2. Pointwise stability and locality. We now discuss our main result, a pointwise 
stability estimate for (2.6) which reveals a generic locality of the TFW interaction. To 
establish this result, we adapt the proof of uniqueness of the TFW equations in [12, 7], 
specialising the class of test functions to 


= 




|V^(a:)| < 7 |^(x)| Vx G 


(3.6) 


for some 7 > 0 . Observe that e G for any 0 < 7 < 7. 


Theorem 3.4. Let mi G and let (mi,0i) Also, let m 2 ■ —)■ M>o satisfy 

II’^ 2 ||l 2 < M' and suppose there exists (m2 ,02) solving (2.6) corresponding to m 2 , 

satisfying U 2 > 0 and 


ll'*^2||rf4^.j(K3) + ||02||i/2^.j(R3) < C{M'). (3.7) 

Further, there exist C = C{M, M', cn), 7 = 'j{M, M', a;) > 0 such that for any f E H,y 

[ ( I'9“'(mi-M 2 )r + 1 ^“'(</’!-<^ 2 )< c [ (mi-ms)^^^. ( 3 . 8 ) 

^|«i|<4 \a2\<2 ^ 

In particular, for any y eM^, 

^ |<9"(mi - U2){y)\‘^ + 1(01 - h){y)\‘^ - ^ f “ iLi 2 )(a:)|^e“^^l*“^l dx. (3.9) 

|a |<2 


Remark 3. Since we do not assume that m 2 E AAl 2 (M',u'), we can not guarantee the 
uniqueness of the corresponding solution (m 2 , 02 )- n 


We can generalise Theorem 3.4 to obtain higher-order pointwise estimates, but this 
requires both inf ui, inf M 2 > 0, hence we need to assume mi, m 2 E M.Hk{M,uj) for some 
k G Nq- 
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Theorem 3.5. Let k eNq and mi, m 2 E Aim {AI, oj). Consider the corresponding ground 
states (mi, 0 i), (m 2 , 02 ) and define 

W = Ui-U 2 , 0 = 01-02, -Rm = 47r(mi - m 2 ). 

Then, there exist C = C{k, M, cu), 7 = 7 (M, cu) > 0 such that for any f E 

[ ( < C' / (3.10) 

^ |ai|<fc+4 |a2|<A:+2 ' \0\<k 

In particular, for any y eM^, 

|R“'w(y)P+ ^ dx. (3.11) 

|oi|<A;+2 |o2|<fc |^|<A; 

Remark 4. It is possible to generalise Theorem 3.4 to treat nuclei described by a non¬ 
negative measure m on satisfying 

sup m{Bi{x)) < M, (HI’) 

xGRS 

and there exist cuq > 0 , cui > 0 such that for all i? > 0 

inf m{Bji{x)) > uqR^ — ui. (H2’) 

The existence and uniqueness of a corresponding ground state {u, 0) is guaranteed by [12, 
Theorem 6.10]. We believe that the arguments used to show [12, Lemma 5.5] and Theo¬ 
rem 3.4 can be adapted to show pointwise estimates similar to (3.8)-(3.9) when mi, m 2 
satisfy (H1’)-(H2’) and that mi—m 2 is absolutely continuous with respect to the Lebesgue 
measure on R^, with a density belonging to 

This result is not sufficient to consider the response of the ground state to a perturbation 
of point nuclei, though it may be possible to treat this using an approximation to the 
identity or by applying similar techniques. □ 

4. Applications 

4.1. Thermodynamic limit estimates. The following result provides an estimate for 
comparing the inhnite ground state with its hnite approximation, over compact sets, thus 
providing explicit rates of convergence for the thermodynamic limit. This is discussed in 
Remark 5. 

Interpreted differently, the result yields estimates on the decay of the perturbation from 
the bulk electronic structure at a domain boundary, generalising the exponential decay 
estimate [ 8 , Theorem 4.6] to arbitrary open hi C R^ and general m E Aii, 2 [M,uj). 

Proposition 4.1. Let m G AiL^{M,u) and {u,(j)) be the corresponding ground state. 
Further, let Q (Z be open and suppose there exists m^ : R^ —)■ R>o such that m^ = m 
on n and ||mf 2|02 < M (e.g., m^ = mxn), then there exists {uQ,(j)n) solving (2.6) 

with m = m^, un > 0 and C = C{M,u)),'y = 7 (M,cu) > 0, independent ofQ, such that 
for all y E Ll 

|5“(« - «p)(2/)l + 1(0 - 0p)(2/)l < (4.1) 

l«l<2 
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Remark 5. Let R> 0 and t C) 0 , then applying Proposition 4.1 with hi = and 

iriQ = gives a rate of convergence for the hnite approximation solving 

(2.4), to the ground state (m,0), 

\\u - URj\w^^^(BR)iO) + \\4> - 4>Rn\\L^{Bn){0) < (4.2) 

This strengthens the result that {uR^,4>Rn) converges to («,</>) pointwise almost every¬ 
where along a subsequence [12]. 

4.2. Locality of the Charge Response. The following result shows that the decay 
properties of the nuclear perturbation are inherited by the response of the ground state. 

Corollary 4.2. Let fc G Nq and mi, m 2 G M.Rk{M,uj). Consider the corresponding 
ground states (mi, 0 i), (m 2 , 02 ) and define 

w = Ui-U 2 , V'= 01 - </> 2 , Rm = 4:7r{mi - m 2 ). 

(1) (Exponential Decay) If Rm £ R^(M^) and spt(Rm) C Br{0), or there exists 7 ' > 

0 such that Yl,\i 5 \<k\^^ ^rn{x)\ < then there exist C = C{k, M,u),'y = 

> 0 depending also on R or 7 ' such that 

|9“iw(a:)| + ^ |(9“^i/>(a:)| < (4.3) 

l“l|<fc+2 \a2\<k 

(2) (Algebraic Decay) If there exist C,r > 0 such that J2\0\<k\^^ < C(l-|-|x|)“'’ 

then there exists C = C{r, k, M, w) > 0 such that 

|5“'w(x)|+ ^ <^^(1 + ^1)-". (4.4) 

|ai|<fc+2 |a2|<fe 

(3) (Global Estimates) If Rm ^ iL^(M^) then there exists C = C{k,M,uj) > 0 such 
that 

II'«^IIh''+ 4(M3) -f ||'0||jyfc+2(K3) < C'||i?m||r/''(K3). (4.5) 

Remark 6. For some of our comparison results, we require only mi G Ai]^ 2 [M,uj) but 
impose weaker assumptions on m 2 . This would not generalise Corollary 4.2 since any of 
the decay assumptions in (1-3) already imply that m 2 G A4j^2{M,uj') for some oj'. 

Remark 7. The estimate (4.3) can be used to study the full non-linear response of the 
ground state to a nuclear impurity. We compare this to the results from the Thomas- 
Fermi (TF) [2, 28, 39] and TFW [18, 35, 40, 29] theories of screening. 

Consider a nuclear arrangement mi G A4i2(M, ca) and model a nuclear impurity at 
the origin with positive charge Z by Zr]{x), where rj G > 0 and J rj = 1. 

Then dehne the perturbed system by m 2 = mi + Zr] G AiL^{Mi,u}i) and consider the 
corresponding TFW ground states {ui, fii) and {u 2 , 02), respectively. From (4.3) of Corol¬ 
lary 4.2 it follows that 

|a“(Mi - U 2 ){x)\ + 1(01 - 02)(a;)| < CZe-^\^\, (4.6) 

l«l<2 

We now compare (4.6) with existing results from the TF and TFW theories of screen¬ 
ing. These models consider the formal linear response (n, V) of the electron density and 
potential to a nuclear impurity at the origin, modelled by the Dirac distribution Zdo, in 
a uniform electron gas. In both models, V satishes the linear equation 

-AV = 47 r[n + Z 60 ], 
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while n solves either the linearised TF or TFW equations. In the TF model, V and n are 
shown to satisfy [28, Page 112], [2, Page 342] 


p-ks\x\ 

V{x) = Z—— 


n{x) 


Zkl 

dvr Ixl 


(4.7) 


where > 0 is a material-dependent constant called the inverse screening length. In the 
TFW model, V and n satisfy [18, 40, 29, 35] 


V(x) = ((a + /3) - (a - , 

Aap\x\ ^ ' 


where a G M, /5 G C satisfy 0 < |/3| < a. The constants a,f3 depend on the material and 
the coefficient Cw, which appears in the dehnition of the TFW energy (2.3). There is a 
critical value of Cy/ below which /3 > 0 and above which f3 is complex, the latter case 
introduces oscillations in the potential and electron density. In either case, both the TF 
and TFW models exhibit screening due to the presence of the exponential term appearing 
in (4.7)-(4.8). 

The lack of a factor of the form in (4.6) can be attributed to using a smeared nuclear 
description for the impurity as opposed to a point description in (4.7)-(4.8). Other than 
this, the similarity of (4.6) to (4.7) suggests that the constant 7 in (4.6) may be interpreted 
as the inverse screening length. In this paper we show there exists 7 > 0 satisfying (4.6), 
however we do not provide any estimates for its value. 

The estimate (4.6) shows that screening occurs in the TFW model, without any approx¬ 
imations made to the model and without any restrictions on the nuclear configurations 
(other than (H1)-(H2)). It should be noted that although (4.6) agrees with existing re¬ 
sults from the TF theory of screening, in metals often the effects of screening are weaker. 
For metals, instead of an exponentially decaying screening factor, Friedel oscillations 
are observed [22, 33, 27]. In this case, the screening factor behaves as |a;|“^/(|a;|), where 
/ : M>o —)■ [—1,1] is an oscillating function and the decay rate r > 0 depends on the Fermi 
surface of the metal. The generic exponential screening factor in (4.6) demonstrates that 
the TFW model significantly overscreens charges. □ 


4.3. Neutrality of defects. An immediate consequence of Corollary 4.2 is the neutrality 
of nuclear perturbations in the TFW equations. This result applies to all nuclear config¬ 
urations belonging to A 4 x, 2 (M,ca). In particular Theorem 4.3(3) strengthens the result of 
[10], which requires mi — m 2 G L^(R^) fl L^(M^) and thus excludes typical point defects; 
see Remark 8 for more details. 


Theorem 4.3. Let mi, m 2 G and define pi 2 := mi — u\ — m 2 -f n 


2 - 


(1) //spt(mi — m 2 ) C BufiO), or there exist C ,7 > 0 such that |(mi — m 2 ){x)\ < 
then pi 2 G L^(M^) and there exist C, 7 > 0 such that, for all R> 0, 




P 12 


< Ce-^^. 


(4.9) 


(2) If there exists C,r > 0 such that |(mi — m 2 ){x)\ < C(1 -|- |a:|) then there exists 
C > 0 such that, for all R > 0, 


'Br{0) 


P 12 


< C(1 + R)2-F 


(4.10) 
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(3) If mi — 772,2 G {e.g., r > 3/2 in (2)) then pi 2 G and 

lim I [ pi2{k)dk = 0, (4.11) 

£^0 |i3e(0)| jBeiO) 

where pi 2 denotes the Fourier transform of pi 2 ■ 

Remark 8 . In a forthcoming article [13], we construct a variational problem to study 
the response of a crystal due to a local defect, using the TFW energy. Arguing as in 
[19], we shall show that any minimising displacement decays away from the defect at the 
rate which corresponds to case (2) with r = 2. In this case (4.10) only provides a 

uniform bound for the charge as opposed to a decay estimate. However, as r > 3/2 the 
global neutrality result (4.11) holds for the relaxed system. 

The neutrality estimates of Theorem 4.3 strengthen those of [10] in the following ways. 
Firstly, our result considers a perturbation of a general nuclear arrangement as opposed 
to a perfect crystal. This allows us, in [13], to consider the response of extended defects 
such as dislocations. In addition, we only require that the nuclear perturbation belongs 
to L^(R^), which we prove rigorously in [13], whereas in [10] the nuclear defect is assumed 
to he in L^(M^) n L^(M^), which fails for typical point defects. □ 

4.4. Energy locality. We now show that the locality result. Theorem 3.5, can be used to 
describe the energy contribution of each individual nucleus. In effect, we will derive a site 
energy potential for the TFW model, which has the surprising consequence that, for the 
study of mechanical response, TFW can be treated as a classical short-ranged interatomic 
potential. Our result gives credence to the construction of interatomic potentials and the 
assumption of strong locality used in hybrid quantum mechanics/molecular mechanics 
(QM/MM) simulations [17]. 

Let p G 0/°(i?j?o(0)) be radially symmetric and satisfy p >0 and p = 1 describe the 
charge density of a single (smeared) nucleus, for some hxed Rq > 0. For any countable 
collection of nuclear coordinates Y = (Yj)j^fq G (M^)^, let the corresponding nuclear 
conhguration be dehned by 

mY{x) = '^^p{x— Yj). (4.12) 

jm 

A natural space of nuclear coordinates, related to the space is 

yL2{M,u) := {Ye (R3)N|my G Ml2{M,u)}. (4.13) 

This space contains many condensed phases, such as crystals containing point defects, 
dislocations and grain boundaries. It does not include arrangements with arbitrarily 
large voids such as surfaces or cracks. However, as the TFW model for surfaces has been 
established [7], it may be possible to obtain locality estimates for surfaces and cracks 
using the TFW model. 

We remark that there exists R' = R'{Ro,uj) > 0 such that for any Y C 3 ^l 2 (M, a;) 

\JBr,{Yj)=R\ (4.14) 

JGN 

For any Y G yL 2 (M,u)) there exists a unique ground state («,</>) corresponding to 
m = my- Naively, we might dehne the energy stored in a region H C R^ by 



(4.15) 
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however, difficulties arise due to the fact that {m — u^) * is not well-dehned. Instead, 
we give two alternative dehnitions for the energy density for an inhnite system: 



Vu ^ + l<Pim - u^). 

(4.16) 

^2{Yr) 


(4.17) 


which both satisfy SiiY ; ■),£ 2 (X !') ^ 

Suppose now that hi C is a charge-neutral volume [44], that is, if n is the unit 
normal to <9^2, then V0 • n = 0 on dVL. Recalling from (2.6b) that 

—A0 = 47r(m — u^) 


we deduce that 

^ f |V0r = ^ f (-A0)0+ f = ^ f (p{m-u'^), 

and hence 


f £i{Y;x)dx= f £ 2 {y',x) dx. 

Jn Jn 

In particular, for hnite neutral systems = m ■ where m G A 4 l 2 (M, co) and 

{ur^, (j)R^) denoting the corresponding solution, the following energies agree on hi = 

^ (^{ttir^ - ulj * (mR„ - 

= - ulj^ 

= + + (4.18) 

We prove this claim is in Remark 12 in Section 6 . Thus, we have derived two energy 
densities, £ 1 , 82 , which are meaningful and well-dehned also for inhnite conhgurations. 

In order to dehne site energies, we require a partition of For each j G N let 
ipj{Y;-) G C^(R^), <yCj(F; •) > 0 satisfying the following conditions: there exist C, 7 > 0 
such that for all Y G yL 2 {M,u) 


'^^j{Y;x) = l, (4.19a) 

jeN 


\iPj{Y-,x)\ and 


dYk 


{Y-x) 




(4.19b) 

(4.19c) 


We propose a canonical construction of such a partition in Remark 9 below. 

Given a family of partition functions satisfying (4.19), we can dehne site energies 


EUY)= f £i{Y-,x)^j{Y-,x)dx, (4.20) 

9k3 

for i = 1,2. A consequence of Theorems 3.4 and 3.5 is that Ej{Y) are local, their 
dependence on the environment of nuclei decays exponentially fast. This is made precise 
in the following theorem. 
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Theorem 4.4. Let i e {1,2}, Y G yL^{M,u)) and {(pj\j G M} satisfy (4.19). Then for 
every k eN, dy^Ej exists and satisfies 


dE^^{Y) 

dYu 




(4.21) 


where C = C{M,u)), 7 = 7 (M, cu) > 0. 


The derivative dy^Ej can be interpreted as the contribution of the atom at to the 
force on the nucleus at Yj. In addition, we show in Section 6.4 that these site energies 
generate the correct total force 


E 


dEj{Y) 

dYk 


_ dE]{Y) 
^ dYk 

jeN 


(fix) 


dmy{x) 

dYk 


dx. 


(4.22) 


Remark 9. Two further canonical requirements on a site energy potential are permuta¬ 
tion and isometry (rotation and translation) invariance. This can be obtained as follows: 

If the partition is permutation invariant that is, for any bijection P : N —?■ N, 

Y o P = (Ypj)j^fq, we have 

fpj{Y o P-^x) = y)p^{Y]x) Vj G N x G M^, (4.23) 

then so are the site energies, 

E]{YoP) = Ei{Y). 

If the partition is isometry invariant that is, for any isometry A : —?■ AY = 

(ylh})jgi!}, we have 

^jiAY- x) = ^j{Y- A-^x) V j G M, X G (4.24) 

then the site energies are also isometry invariant, 

E}{AY) = E]{Y). 

Both statements are proven in Lemma 6.9. 

A canonical class of partitions satisfying (4.19) as well as (4.23), (4.24) can be con¬ 
structed as follows: Let tp G C^(M^), ^ > 0, be radially symmetric and satisfy 

|^(x)| + |V^(x)| < (4.25) 

(p{x) > c > 0 on Pjj/+i(0). (4.26) 


For example, this holds for ^(x) = and for standard mollihers with sufficiently 

wide support. 

Then, for j G M, we can dehne 


y)j{Y]x) 


Tjx - Yj) 
XlyeN ~ ^j') 


(4.27) 


It is easy to see that this class of functions are well-dehned and satishes all requirements. 

□ 


Remark 10. Alternative constructions of energy partitions include Bader volumes and 
charge-neutral volumes [4, 44, 32]. Bader volumes partition space into regions such that 
the flux of the electron density on the boundary is zero, while charge-neutral volumes are 
dehned so that each region has zero charge. The construction of these volumes is not 
unique, like our dehnition of a partition. Bader volumes were developed as a means to 
dehne atoms within molecules [4]. 
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With this in mind, using a partition we may assign a portion of the electron density 
to each nucleus in the system. We refer to a nucleus paired with it’s associated partition 
of the electron density as an effective atom. Due to the screening that occurs in the 
TFW model, the interaction of two effective atoms decays exponentially as the distance 
between the nuclei grows. In comparison, the interaction of two neutrals atoms separated 
by a sufficiently large distance r in the TF model has been shown to decay at the rate 
r~^ [9]. This suggests that due to the overscreening of the TFW model, the interaction 
of the effective atoms is considerably weaker than is realistic. However, for the purpose 
of simulating quantum systems, in particular applying the strong locality principle [17], 
the weak long-range interaction of the TFW model is a desirable property. □ 

Remark 11. The estimate shown in Theorem 4.4 is a theoretical result which can be 
applied to simulate defective crystals, though we do not pursue this. Locality estimates 
have been established for the tight-binding model and subsequently used to construct 
QM/MM hybrid methods [14, 15]. □ 


5. Conclusion and Outlook 

The two main results of this work. Theorems 3.4 and 3.5, are stability and exponen¬ 
tial locality estimates for the TFW model, which apply to general condensed nuclear 
conhgurations. 

We have demonstrated in Section 4 that it can be used to extend and strengthen a range 
of existing results on the TFW model. A particular strength of our results is that they 
apply to general nuclear configuration in Ai u), whereas the previous analyses of the 
TFW model have focused on (near-)crystalline arrangements or the homogeneous electron 
gas. This generality will be valuable when exploring the consequences of our analysis for 
studying models for the mechanical response problem in [13], where we generalise [19] to 
electronic structure models. 

A further application, that we will develop in a forthcoming work is a study of the 
Yukawa potential as a model approximation [36]. Adapting Theorems 3.4 and 3.5 we can 
consider the difference between the Coulomb and Yukawa ground states for a given nuclear 
conhguration and prove uniform error estimates in terms of the screening parameter in 
the Yukawa model. 

Two key difficulties in the analysis of electronic structure models are (i) the exchange 
and correlation of electrons due to the antisymmetry of the electronic wavefunction; and 
(ii) the interaction of charged particles (positive nuclei and negative electrons) via the long- 
range Coulomb potential. Since the TFW model is orbital-free it does not account for (i), 
however it fully incorporates Coulomb interaction. In this regard it is perhaps surprising 
that the TFW model satishes the extremely generic locality property we obtained in 
Theorems 3.4 and 3.5. 

The Hartree-Fock and Kohn-Sham models take both effects into account and whether 
these models permit a similarly strong notion of locality is an open problem. It is clear, 
however, that such results cannot be obtained in the generality that we considered in the 
present paper. Since charged defects exist in the reduced Hartree-Fock model [11] and as 
locality implies neutrality, this suggests that a locality property cannot hold for general 
condensed phase arrangements in the reduced Hartree-Fock model, which is the simplest 
model in the Hartree-Fock/Kohn-Sham class. 


LOCALITY OF THE THOMAS-FERMI~VON WEIZSACKER EQUATIONS 


13 


6 . Proofs 


This section contains the proofs of the main results. Proofs of results in Sections 3.1, 
3.2 and 4 are found in Sections 6.1, 6.2 and 6.3 respectively. 

The following is a preliminary result used in the construction of the space 

Lemma 6.1. Suppose m : —)■ M>o and m G then (H2) is equivalent to the 

following statement: there exist ojq^uji > 0 such that for all R > 0 

inf / m{z) dz > cuqR^ — uii- (6.1) 

JBaix) 

Proof of Lemma 6.1. Clearly, (6.1) implies (H2), so suppose m satisfies (H2), then there 
exists > 0 such that 

inf / m{z) dz > 1. 

For R > 0 and x' G M^, let Qr{x') C denote the cube of side length 2R centred at x\ 
which contains Br{x'). Also, let i?2 = V^Ri, which ensures that Br^(x) D Qrj^{x) for all 
X G Further, let R> R 2 , hence R = kR 2 , for some k>l. Then 

inf / m{z) dz = inf / m[z) dz > inf / m{z) dz 

Jbr{x) XG1R3 

> \_k\^ inf [ m{z) dz > [k\^ 

^'GR3 

/kV 

K2) = 875 =^"“''' 

Now define Ui := ojqRI > 0, then it follows from (6.2) that (6.1) holds for all R > 0. □ 


6.1. Proofs of Uniform Regularity Estimates. The following lemma features in the 
proofs of both the existence and uniqueness of the TFW equations and is found in [12]. 

Lemma 6.2. Let a G fl L°°(M^), then define the elliptic operator L = —A + a. 

Suppose that there exists u G satisfying u > 0 and Lu = 0 in distribution. Then, 

the operator L is non-negative, that is for all (p G R^(M^) 

{(p, Lp) > 0. (6.3) 


The proof is shown in [12] but is included here for completeness. 

Proof of Lemma 6.2. Let R > d and dehne fl = i?i?(0) and consider L as an operator 
on L^(f2) with domain H‘^{Q) fl Hq{Q). Then L is a self-adjoint operator with compact 
resolvent hence has a purely discrete spectrum. Since a G iLiQj,(R^) it follows that the 
smallest eigenvalue Ai(f2) is simple and has a positive eigenfunction vn G Hoi^) [23, 
Theorem 8.38]. In addition, by standard elliptic regularity vq G ^ [ 20 ] 

and solves 


(-A + a)v^ = Xiipdjv^. 


Testing this equation with u and using integration by parts, we obtain 


(6.4) 


As uq > 0 on fl and vq_ vanishes over clfl, it follows that <0. It follows that the 
left-hand side of (6.4) is non-negative, hence Ai(fl) > 0 . As this holds for fl = BR{d), for 
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any i? > 0, we deduce that for all ip G {ip, Lip) > 0. Using that a G and 

the density of U^(R^) in ih^(R^), it follows that for all ip G {ip, Lip) >0. □ 

We now show uniform estimates for hnite systems corresponding to truncated nuclear 
distributions. This result is essentially [12, Proposition 3.5], however as we require uniform 
regularity estimates, we provide a complete proof. 


Proposition 6.3. Let m : R —)■ R>o satisfy 

||iH||i2^.^(K3) < M, (6.5) 

and Rn t oo, then define the truncated nuclear distribution = rn-XBu^io)- The unique 
solution to the minimisation problem 


rTFW/ 


(mR„) = inf { 


V G > 0, 


v^ = 


mRu 


yields a unique solution to (2.4) 

- (t)R„UR^ = 0 , 

= 47r(mR^ - w^J. 

which satisfy the following estimates, with constant C independent of Rn: 

+ (67) 

ll'(’Rnllw^„i((E3) 7 0(1 + (6-8) 

Proof of Proposition 6.3. If m = 0, then for all Rn, clearly ur^ = (pR^ = nriR^ = 0 satishes 
(2.4) and (6.7)-(6.8). 

If m ^ 0, then there exists a constant ^ 0 such that Rn > Ro ensures that 
/r 3 mRn > 0. Recall 

= J I Vn|^ + J + ^D{mR„ - mR^ - n^). 

For each Rn, consider the minimisation problem 


rTFW/ 


(m^J = inf U™(n,mRj 


V G H^(R^),v > 0, 


mRn 



The constraint = f ^3 mR^ ensures the system is charge neutral, and by [31, Theorem 

7.19] there exists a unique non-negative minimiser ur„ G H^{R^) to {mR„) solving 

- (("^Rn - ur„ = -Or^ur^, (6.9) 

/ “L = / > 0- (6-10) 

JK3 JR3 

Here 9r^ > 0 is the Lagrange multiplier associated with the charge constraint (6.10) 
[31, 12]. Dehne (pR^ : R^ —)■ R by 

0Rn = {{mRn - UrJ * h) “ ^Rnl 
so we can express (6.9) as the Schrodinger-Poisson system (2.4) 

-^URn + ^Ur! - (t^RnURn = 0 , 

= 47r(mR,, -RrJ- 


( 6 . 11 ) 
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Decompose 

{rriR^ = {niR^ - * (^Xbi(o)) + “ “L) * (n^Bi(o)=) > 

then as ur^ G ^ and m e applying Young’s inequality gives 


{rriR^ -ulj*n 


< Wil^Rn - UrJ\\l5/3(R3) 


llXBiiO) 


L5/2(R3) 


II - “iJ„)llL 7 / 5 (R 3 ) 


-\XBrioy 


LY2(R3) 


< C + -Rn^^'^)ll’^i?nlU2(R3) + ||MR„||^1(R3) 

^ C (i^n^ + -^n^'^'^)ll’^llLLif('*h + ll“lJnllHl(R3) 


<C ((«=''=+ fi«OAf+||«„, 


I Hi 


By [30, Lemma 11.25] we deduce that {rriR^ — u\^) * is a continuous function vanishing 
at inhnity. It follows that (pR^ G and is also continuous. Also, |V0 _r„| G 


1 

Svr 


\^4>rI = 


Stt 

\l 


4>Rn (-^0Rn) 
83 


^ / {ruR^ - ujij 
^ JR3 

{4>r„+0rJ {mR„ -ulj 

-ulj* {niR^ -ulj, 


( 6 . 12 ) 


hence (j)R^ G Now, consider ur^ G which solves 

-Aur^ = ~u£ + (t>R^URr.- (6.13) 

The right-hand side can be estimated in L^(M^) by 

K 1 /“^ 5 7 /T 

lls^Hn “ 0Hn'WHnlU2(R3) < -\\Ur^ I|l 2(R3) + | [ 11 L^o (R3 ) 11 11 ^2 (r3 ) 

< C'||-Ui?„||^f(jj3) + |l(/)iJ^||ioc(R3)||Mi?„||/^l(R3), 

which implies ur„ G iL^(R^) as Aur^ G By the Sobolev Embedding Theorem 

[20] ur^ G hence ur^ is continuous. Also, by [5, Lemma 9], 

ur^ decays at inhnity. We now justify this. Recall (6.13) and since ur^ > 0, we have 
-Aur^ < (I)r„ur„, hence 


-Aur^ -h ur^ < (1 + 0h„)mh„ 

As (j)R„ G L°°(M^) and ur„ G iL^(M^), the right-hand side belongs to L^(M^) hence by the 
Lax-Milgram theorem there exists a unique qr^ G satisfying 

-Aqr^ + Qr^ = (1 + 0 h„)wr„- 

Moreover, using the Green’s function qr^ = * (1 - 1 - (t>R^)uR^ and since (1 -|- 

4>r„)ur^ G L^(M^) by [30, Lemma 11.25] gR^ is continuous function that decays at inhnity. 
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hence G It follows from the comparison principle that ur„ < qr^, so 

ur^ G and decays at inhnity. 

Using that UR^,(j)R^ + 9r^ are continuous and decay at inhnity, by arguing as in [41], 
there exists a universal constant Cs > 0, independent of the nuclear distribution, satisfy¬ 
ing 

0 < < Cs, (6.14) 

+ Cs. (6.15) 

As ur^ > 0, from the Solovej estimate (6.15) we obtain the uniform lower bound 

> -Cs. (6.16) 

We aim to show a uniform upper bound for 4>r„, which together with (6.15) will yield the 
uniform estimate 

+ ||0r„||l°°(k3) < C{M), (6.17) 

which is independent of Rn- 

If (f)R„ is non-positive, then (6.17) holds as 

+ ||0R„|U°°(R3) < 2 U 5 . 

Instead, suppose that (j)~^ is non-zero at some point in M^. By (6.11) (pR^ is a continuous 
function that converges to a negative limit at inhnity, G C'c(M^), hence there exists a 
point Xr^ G such that 

(pR„i^Rn) = II0 rJU“(RU > 0 . (6.18) 

Without loss of generality, we assume = 0. 

We now show that ur^ > 0 on M^, arguing by contradiction. Suppose that there exists 
z G such that ur^{z) = 0. Since ur^ is a non-negative, continuous function decaying 
at inhnity, there exists i/n G such that 

URAVn) = sup ur„{x). 

3;G1R3 

Let R > lun — z\, then by the Harnack inequality [42], we infer 

0<UR„{yn)= sup ur„{x)<C{R) inf ur^{x) = ur^{z) = 0, 

x£Bji{yn) 

so ur„ = 0. This contradicts the charge constraint (6.10) /r 3 = /r 3 > 0, hence 

uRu > 0 on 

As ur^ G H^{R^)r\L°°{R^),4>R^ G H Cf {R^)r\L°°{R^) and ur^ > 0, Lemma 6.2 implies 
that Lr^ = —A -|- — (j)R^ is a non-negative operator. 

Choose ip G C“(i?i(0)) satisfying 0 <<y 5 <l, 9 ? = lon i?i/ 2 ( 0 ), /r 3 9 ?^ = 1 and 
/r 3 =: c^, then for y G dehne ipy G C^{Bi{y)) by ipy = ip{- - y). As Lr^ is 

non-negative (6.3) implies 

Lr^^v) = [ iVifyl"^ + / ( 

jR3 JR3 V-J / 

which can be re-arranged and expressed using convolutions as 

I (“K* * ^ (^^Rn * “ ^3 

= (0R„ - W) + 


(6.19) 
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Observe that (j)R^ * 99^ solves 

= 477 {niR^ * 

We estimate the hrst term using (6.5) 


* 




( 6 . 20 ) 


(mR„ * Lf^) {x) = / mR^{y)Lp^{x - y) dy 

JBi{x) 

< / m{y) dy < Oo||m||i2 ^s) < CqM. 
J Bi{x) 


( 6 . 21 ) 


For the second term, using the convexity of t i-A for f > 0 and the fact that J = 1, 
applying Jensen’s inequality and (6.19) we deduce 

dTTul^ * if^{x) > * (p2(x) 


> 


- y)i/{y) dy 


u 


,4/3^ 

A. 


\ 3/2 

(x-y)j ip^(y) dy 


\ 3/2 

dy] 


( 6 . 22 ) 


= ^ i^Rr^ W) + ^ • 

Combining the estimates (6.20) - (6.22) we conclude that 

(0R„ * <d^) + (0Rn * < CqM. 

By (6.11), as (j)R^ is a continuous function that converges to a negative limit at inhnity, 
(pR^ * 99^ also shares these properties. Now consider the set 

S' = { X e I > 0 }, 

it follows that S is open and bounded and that pR^ * 99^ — = 0 on dS. Observe that 

the constant function h = {CqMY/^ satishes 

- Ah + = CqM on 

0 = < h in dS, 

so by the maximum principle pR^ * < (c^, + over S', and also on S'^, hence 

0R. *9^'<C'i(1 + M2/3), (6.23) 

where Ci = max{c^, is independent of M. 

Observe that 

* ‘d' = ^R^ * + 0«„ *y^^<Cs + C,{l + M2/3) = 0(1 + 

and by estimating (2.4b) directly, that 

-^01 = -^0AnX{</«„>o} = 477 {mR^ - X{0«„>o} < 477mK„X{0«„>o} < 477mj 

As 0 < <y9 < 1 and 99 = 1 on i?i/2(0), then 


^Rr. 


^ Bi/2{0) 


pR„{x) dx < {p-^^ * ip^) (0) < 0(1 + M^/^). 


(6.24) 
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Using a change of variables, (6.24) can be expressed as 

.1/2 


/ dx = 

Define / : [0,1/2] —)■ M by 

m = 

and snppose that for all t G (1/4,1/2) 

/(t)>4^ 

then 


/o JStiO) 


d^t(7) dt. 


St{0) 




1^1/2(0) 


dx > 0, 


' Bi/2{0) 


(I>rAA dx = 


.1/2 


.1/2 


fit) dt> fit) dt > 

7i/4 


^^1/2(0) 


dx. 


which gives a contradiction, hence for some t G (1/4,1/2) 


JSt{0) 

Since t > 1/4, (6.25) implies 


^L(*7) dS,(7) < 4 / .PiJx) di < C(1 + 


/i3i/2(0) 


/ d^i(7) = [ 01(^7) d^t(7) 

Jst{o) PnUJI Jst{o) 


< 


l^i(0)| 

C + M2/3 


SdO) 

< Cil + M^^A =■ CiiM). 


\Si/A0)\ 

We now construct an upper bound for as follows. Let (fi satisfy 

-A01 = 0 in BtiO), 

As (pi is harmonic, it satishes the mean value property 
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(6.25) 


0i(O) < f 0R„(t7) d^t(7) < CiiM). (6.26) 

JSt{ 0 ) 

Then consider the Dirichlet problem 

—A02 = 47rm in BtiO), 

(p 2 = 0 on S't(O). 

By Lax-Milgram, this has a unique weak solution 02 G H^iBiiO)). By standard elliptic 
regularity theory [20] 02 € iL^(i?0O)) ^ iBtiO)) and 


1102 


^2|Ic0.i/2(s70)) - 
The constructed functions 0i,02 satisfy 

-A0;^^ < -A(0i + 02) in BtiO), 


C||02||r2(5j(o)) < U||m||i2(5^(o)) < Ut^^^||m|]^2^,^(]R3) < CM. (6.27) 


Ar„ - 01 + 


on 


StiO), 


hence by the maximum principle 0]^^ < 0i + 02, in particular (6.26)-(6.27) imply 
I|0 /J„||l°°(K 3) = 0ij„(O) < 0l(O) + 02(0) < U(1 + M), 
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where the right-hand side is independent of Rn- Combining this with the lower bound 
(6.16) and the Solovej estimate (6.15), we obtain the estimate (6.17) 

+ ||0R,i||l°°(E3) < <^(1 M). 

It follows immediately that for all a; G and r G [1, cxo] 

\\uR„\\Lr{B2{x)) < C{1 + (6.28) 


independently of both x, r and Rn- Using (6.17) and (6.28), we now obtain uniform local 
estimates for the right-hand side of (6.13) 

= —Ujll + (t^RnURn 
by 

- 4>RnURAL^{B2{x)) < C\\^U^^ - 

< + ||0A„||l“(]R3)||-Ur„||i,oo(r3)) 

< C(1 + M^/^). 


Consequently, for any x G M^, the elliptic regularity estimate [20] gives 

\\urJ\H'^(Bi{x)) < “ 0R„'«^A„IU2(B2(a;)) \\u rJ\ (B2{x))) 

< C(1 + M^/^) + C(1 + ^ 29) 


As (6.29) is independent of x G we obtain 

ll'^A„llH2^if(R3) < C(1 + (6.30) 

Applying a similar argument to estimate the right-hand side of (2.4b) 

= 47r(mij„ - 


yields (6.8) 

ll0Rn|lH2^ij(R3) < <^(1 + 

Using that (j)R^ G and arguing as in (6.30), we obtain the desired estimate (6.7) 


ll'^An llH4nif(R3) < C{1 + 


□ 


We remark that while the constants appearing in the hnal estimates (6.7)~(6.8) depend 
on c^, they are independent of M. 


Remark 12. We now justify the claim that for hnite and neutral systems and for 11 = 
the three energies shown in (4.18) agree. Recall (6.12), which shows that the Coulomb 
energy can be expressed as 

\ {i^Rn - mL) * h) - “D = ^ “ ^rJ ^ ^ ^3 I I ^ ’ 

so it follows that the energies dehned in (4.18) agree for fl = M^. □ 


We now discuss passing to the limit in (2.4) to obtain regularity for the inhnite system. 

Proof of Proposition 3.1. First suppose that spt(m) is bounded, then for sufficiently large 
Rn, m = mjR and hence by Proposition 6.3 (n, </>) = {uRnif’Rn) solves (2.6) and satishes 
the desired estimates (3.2)-(3.3). 

Now suppose spt(m) is unbounded, then the estimates (6.7)-(6.8) of Proposition 6.3 
guarantee that the sequences ur^, (j)R^ are bounded uniformly in Consequently, 
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there exist m, 0 G fl such that along a subsequence UR^,(j)R„ converges 

to u,(j), weakly in strongly in for all i? > 0 and pointwise almost 

everywhere. It follows from the pointwise convergence that n > 0 and 

||ii||i,°°(r3) < C {1 + 

||0||l°°(r3) < C {1 + M ). 

Passing to the limit of the equations (2.4) in distribution, we hnd that the limit (u, 0) 
solves 

5 

—Am + — 0M = 0, 

o 

—A(/) = 47r(m — u^). 

Arguing as in (6.7)~(6.8), we deduce that the desired estimates (3.2)-(3.3) hold 

||0||J‘^(r3)<C^(i + m3/2). □ 

Proof of Proposition 3.2. As m G M.l 2 {M,uj), it satishes (H1)-(H2), hence by Theo¬ 
rem 2.1, the solution (w, 0) of (2.6) dehned in Proposition 3.1 is unique and satishes 
inf M > 0. Now suppose 

inf inf u(x) = 0, (6.31) 

mGAI^2 xGR® 

we show that this contradicts the assumption that for all m G M.l 2 {M,uj) and R> 0 


inf 

3:eR3 


Bnix) 


m[z) dz > UqR^ — oji. 


It follows from (6.31) that there exists rrin G A4 l2 (M, ca) with corresponding solution 
(m„, (fn) and Xn G such that for all n G N 

/ N 1 

'^nipjn) A • 
n 

Recall the uniform estimates (6.7)-(6.8) from Proposition 3.1 

ll'nn||ir4^ij(R3) < C{1 + 
ll</’n||ir2^ij(R3) < C(1 -I- M^'^^). 

It follows that 


5,, 4/3 
3 


Un 


^unif' 


< c\\u 


14/3 


+ \\4>n 


’*IIL°°(R3) ' llYn||L2^.j(R3)|lMn||L°° 


< C{M). 


As \un'^ - 4>nUn £ G i7,)jjif(M^) and «„ > 0 solves 


LnUn ■= I -A - 


Un 0 , 


applying the Harnack inequality [42], and observing that the coefficients of are uni¬ 
formly estimated by Proposition 3.1, this yields a uniform Harnack constant, hence for 
all R > 0, there exists C = C{R, M) > 0 such that 

C 

sup Un{x) < C inf Un{x) < —. 
xeBnixn) xeBnixr,) n 


It follows that the sequence of functions Un{- + Xn) converges uniformly to zero on compact 
sets. Consider the ground state (m„, 0„) corresponding to the nuclear distribution m„. 
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Recall that solves the following equation in distribution 

-A(j)n = dvr (m„ - ul) . (6.32) 

We translate the system and then pass to the limit in (6.32) as n tends to inhnity. To do 
this, we use the following estimates, which are translation invariant: 

||"ln(- +a:n)||L2^.j(K3) < M, 

||0n(- + ^n)||r/2^.j(R3) < C{M). 

It follows that, up to a subsequence, 0„(- + x„) converges to <p, weakly in 
strongly in for all i? > 0 and pointwise almost everywhere. Moreover, m„(- + 

Xn) converges to m, weakly in L^(i?j^(0)) for all i? > 0. By applying the Lebesgue- 
Besicovitch Differentiation Theorem [21] we deduce that rh G Passing to the 

limit in 

-A0„(- + Xn) = 471 (mn(- + Xn) - M^(- + X„)) , 
it follows that 0 is a distributional solution of 

—A0 = 47rm. (6.33) 

Arguing as in [12, Theorem 6.10], we show that for all i? > 0 


/ m{z) dz < CR. 

Jb aio) 

As m G A4i2(M,cu), this leads to the contradiction that for all i? > 0 


(6.34) 


OJqR — UJ\ ^ 


'Hh(O) 


m{z) dz < CR. 


To show (6.34) choose G C“(R2(0)) such that 0 < (p < 1 and (p = 1 on i?i(0). Let 
R > 0, then testing (6.33) with (p{-/R) gives 


I (l){z){A^){z/R) dz = in 


' B2r{0) 

The left-hand side can be estimated by 
1 


' B2r{0) 


m{z)ip{z/R) dz. 


(6.35) 


R2 


' B2r{0) 


(j){z){Aip){z/R) dz 




|-^2ij(0)| 

R2 


< CR, 


where the constant C* > 0 is independent of R. As m > 0, from (6.35) we obtain (6.34) 



m{z) dz < 


' B2r{0) 


m{z)ip{z/R) dz < CR. 


The contradiction ensures that there exists a constant cm,uj > 0 such that for all m G 
A4i2(M, cn), the corresponding electron density u satishes 


inf u(x) > Cm oj > d. □ 

a;GR3 


Proof of Corollary 3.3. Our aim is to show by induction that for all k G No, if m G 
then the corresponding solution (u, 0) to (2.6) satishes 
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In Proposition 3.1, by combining the estimates (3.2) and (3.3), it follows that (6.36) 
holds for the case /c = 0: for all m G AiL^{M,u)) the corresponding solution {u, (j)) satishes 

+ ll0llH2^if(lR3) < C{M,Uj). 

We now show the induction step. Suppose the result is true for k G No, then consider 
m G A4Hk+i{M,u)) C so by the induction hypothesis the corresponding 

solution (m, (j)) satishes 

ll'^llg,(^ + hR3) + II‘/’II_h'^+2(]r3) ^ C • (6.37) 

We remark that as 0 < cm,w C{M) and u G it follows that for all r G M, 

G i7^+f(R3). As (m, 0) solve (2.6) 

— An = + 4>u, 

o 

—A(/) = 47r(m — n^), 
by standard elliptic regularity theory [20] for any x G 

I|0||h'=+3(Bi(x)) < C (||m — U^\\h>=+^{B2(x)) + II0 |It2 (B2 (x))) 

< C {\\m\\Hk+i(B 2 {x)) + \\u^\\Hk+i(B 2 {x)) + |I0|U2(B2 (x))) 

< c (|I"^IIh^+)(r3 ) + ll0lL2^if(R3)) +c (^k + 1 , ||n||^fc+j(],j3)) 

< C||m||^fe+j^]g3^ + C {k + 1, ||m||jyfc^.^(]K3),cjj 
<C {k + l, ||m||^fc+j(]g3),a;j , 

hence 

||0lln^lj(R3) “ ll‘/’lliL*^+3(_Bi(x)) A C + 1, ||?TT.||. (6.38) 


We use an identical argument together and apply the estimate (6.38) to deduce 



Combining (6.38) and (6.39) we obtain the desired estimate 

II'*^IIh*+®(r 3) + ll'!^llir^;t3(M3) < C {k -\- , 

which completes the proof of (6.36) by induction. 


(6.39) 


□ 


6.2. Proofs of Pointwise Stability Estimates. To prove Theorems 3.4 and 3.5, we 
adapt the proof of uniqueness of the TFW equations, shown in [12, 7]. Due to the length 
of the argument, we shall prove several intermediate results. Before showing these results, 
we outline the structure of the proof. 

First, we state two alternative sets of assumptions on nuclear distributions mi, m 2 '. 
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(A) Let fc = 0, mi G and let denote the corresponding ground 

state. Also, let m2 : —)■ M>o satisfy ||m2||L2 < M' and suppose there exists 

(m 2,02) solving (2.6) corresponding to m2, satisfying M2 > 0 and 

+ 110211^2 (IR3) < C{M'). (6.40) 


“2 

unit 

In addition, we assume that either m2 ^ 0 and U 2 > 0, or m2 = M2 = 02 = 0. 

(B) Let k G No, mi,m 2 G and let (mi,0i), (m2,02) denote the correspond¬ 

ing ground states. (Note that (B) implies (A), with M' = C{M).) 

Remark 13. We point out that in (A) we assume M 2 > 0, while in Theorem 3.4 we only 
require M2 > 0. The restriction M2 > 0 allows us to directly use results from [12], in 
particular Lemma 6.2, and will be lifted via a thermodynamic limit argument in the third 
part of its proof on page 32. 

Throughout the remainder of the paper we use the notation 

w = Ui - U2, -0 = 01-02, Rm = 47 r(mi - m2). 

By treating the coupled system of equations as a linear system and by exploiting the 
coupling between the electron density and electrostatic potential arising from the Coulomb 
energy term of the TFW functional, we obtain the following initial estimates 

Lemma 6.4. Suppose (A) holds, then there exists C = C{M,M',u) > 0 such that for 
any f G 


{w^ + \Vw\^+\V'l/j\^)e<C( Rm'll^e+ (lT2 + 02)|Ve|‘ 


(6.41) 


To control the 0-dependence on the right-hand side of (6.41), we require an estimate 
of the form 


re < c { / / { w ^+ r ) m - 


(6.42) 


which holds for ^ G Hi, i.e. f G iL^(M^) satisfying |V^| < ^ on M^. 

Suppose (6.42) holds, then applying Holder’s inequality and (6.41) yields 


/ < C ( / Rie+ / (M;^ + 00|Ve| 

JR3 VJR3 Jr3 

To remove the term J(w‘^ + 0^)|V.^P on the right-hand side, we simply restrict from 
^ G iLi to a narrower class of test functions, 

H, = {^e I \Vf{x)\ < 7 |e(a;)| Vx G M^, 

where 7 = min{l, (2C")“^/^} > 0, to show 


(m;2 + I Vm;^ + 0^ + IV0ne^ < 2C" / Rie 


(6.43) 


In order to show (6.42), we adapt the argument used in [7] . At the same time, since the 
equations for (^,0) hold pointwise, we obtain additional estimates for Ate, A0. 

Lemma 6.5. Suppose (A) holds, then there exists C = C{M, M',u),'y = 'y{M, M',u) > 0 
such that for any f G 


+ |Vm;|^ |Am;|^ 0^ -t- |V0|^ + |A0|^ < C 


Rie- 


(6.44) 
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Clearly Lemmas 6.4 and 6.5 hold also under the assumption (B) since (B) implies (A), 
with M' = C{M). In the case (B) where mi, m 2 are both uniformly bounded below and 
have higher regularity, arguing as in Corollary 3.3 and Lemma 6.5, we obtain improved 
estimates for w and ^|J. 

Observe that in Case (B), M' = C{M). Due to this, we omit the dependence of M' in 
the constants that appear in the following lemmas, whenever we assume (B) holds. 

Lemma 6.6. Suppose that either (A) or (B) holds, then there exist C = Ca{M, M', cn), 7 = 
'jAiM, M',u) > 0 or C = Csik, M,uj), 7 = > 0, where 'Jb independent of k, 

such that for any f G 

[ ( (6.45) 

^ |ai|<fc+4 \a2\<k+2 ' |^|<fc 

In particular, for any y 

\d°'^w{y)\‘^ + ^ < C f dx. (6.46) 

|ai|<fc+2 |a2|<fc \y\<k 

We remark that in the following proofs, all integrals are taken over unless stated 
otherwise. 


Proof of Lemma 6.4- Case 1. First suppose that m2 ^ 0 and U 2 > 0. Recall that mi G 
hence by Propositions 3.1, 3.2 and (6.40) 

l|w2|k4^.j(R3) + I|02 ||h2^.^(R3) < C(M'), 
inf uAx) > Cm ui > d. 

By the Sobolev embedding: for all k E Nq and x eM.^ H^^‘^{Bi{x)) C'^’^/^(i?i(a;)), so 

it follows that 


||Mi||iy2.°°(M3) + ||0 i||l°o(M3) < C{M), 

||m2||iT2.°°(R3) + ||02||l°°(R3) < C{M'), 

hence w = U 1 —U 2 E iLy^;f(M^) O hF^’°°(M^), fj = (pi — <p 2 E H Cf {R^) nL°°(M^), and solve 


-Aw = ^ j + (piUi - (P 2 U 2 , 

-Alp = AtT [ul - Ml) + Rni, 

pointwise. Let f E then test (6.47a) with to obtain 

j Vw ■ V{wf^) + ^ y (“ 1 ^^ - - j {(plUl - (p 2 U 2 )wf‘^ = 0. 

We will use the following rearrangements 


(6.47a) 

(6.47b) 

(6.48) 


0iMi - (P 2 U 2 = (6.49) 

j Vw • V«2) = j |V«)|2 - j (6.50) 

J • v{f^e) = j\^m\" - J ( 6 . 51 ) 
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To estimate the second term of (6.48), observe that Proposition 3.2 and (A) imply that 
inf Ml > cm,u) > 0 and recall the assnmption -U 2 > 0. It follows that for f + 

> 0 

— M2 ^^)(Mi — U 2 ) = {u\^^ + U2^)w^ + 

> (Mi^^ + U2^‘^)w‘^ 

> + U2^)w'^ + uw"^. (6.52) 

Combining the estimates (6.48)-(6.50) and (6.52), we obtain 

I |V(wOP + ^ J + U2^^)w‘^^'^ - ^ + ^ j 

< j j ^p{ul - ul)i‘^. (6.53) 

We dehne the following operators 

Li = -A + - 01, 

L 2 = -A + - 02, 

L = -Li + -L 2 = —A + ^ 

As Ml, M 2 > 0, Lemma 6.2 implies that ^ 1,^2 are non-negative operators, hence for any 

e 

{ip, Lp>) = ^((^, Lip)) + ^((^, L2P>) > 0. (6.54) 

Observe that as m; G 1T^’°°(M^) and ^ G iL^(M^), G iL^(M^). We can express (6.53) as 

(m;^,L(m;0 ) + J < J + ^ J ^p{ul-ul)^^. (6.55) 

To control the hnal term of (6.55), we begin by testing (6.47b) with 0,^^ to obtain 

J V0 • V(0e') = 471 y 0 (m2 - uDe + J Rmi^e- (6.56) 

Rearranging (6.56) and applying (6.51) yields 

i j ij(ui - =A y _ i_ y 

= A / _ i_ y |V(«)|'^ + A y (6.57) 

Combining (6.55) and (6.57) yields 

{w^,L{w^)) + u J + ^ j |V(0OP 

< ^ / ^/0^|ve 


(6.58) 
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As = V('0O ~ V'V.C, we have 

[\vi:\^e<c[ I iv(v^or+ / 


<ci^j Rmi^e+jiw^+rmi" 

Combining the estimates (6.58)-(6.59), we obtain 

K,LK))+i. j w^e+j 

<c( [ R^^pe+ 


(6.59) 


(6.60) 


Next we obtain an estimate for J | using the fact that L is a non-negative operator. 

We can express L as 

(j)l + 4>2 tt2 /in,S\ 

L = -A a, where a = -----G 

6 / 

From (6.54), we have shown that L = —A -|- a > 0 in the sense that (</?, L(f) > 0 for every 
ip G So for e G (0,1) 

A = (1 — £)(—A -|- o) -|- £(—A) -p ECL ^ £(—A) — sjlcill^o 

Applying this to (6.60) gives 


^ / |V(wOP + {i^- £\\a\\L- 


2c2 


<C[^J R„i-e+ + 

so choosing e = min{ ^||^|j'^^ , ^}, we deduce 

J \V{wO\^<c(^J Rm^pe + + 

and since .^Vtc = V{w^) — wV^, we deduce 

J\vw\^e<cQ Rmi^eR (o.oi) 

We combine the estimates (6.60) and (6.61) to obtain the desired estimate (6.41) 

fw^e+ f |Vu;pe^+ [\Vm^<c( [ Rmi^e+ + 


and observe that this estimate is valid for any ^ G 

Case 2. Suppose now that m 2 = M 2 = 02 = 0, then the argument used to show (6.53) 
holds to give 


|V(wOP+^ J J + J 

< f + ^ f 


(6.62) 
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Now using that Li is a non-negative operator, we obtain 

^ j |V«)|^ + z/ J < ^{(p,Li(p) + ^J |V«)|^ + i^ J 

= J |V(wOP + ^ j J (piw^^^ + u J 

< J + ^ j 

Then applying the estimates (6.56)-(6.60) yields the desired estimate (6.41): for all ^ G 
H\R^) 

+ j j\Vm"<c(^j Rm^e + + □ 


Proof of Lemma 6.5. To obtain an integral estimate for ip, hrst recall (6.47a), that w 
solves 

A , 5 / 7/3 7/3 \ 01 + 02 Ui + U2 . 

-Uf J - - - W= - - -0, 

then testing this equation with 0^^, for ^ G yields 


III + / 2^2 [a / ^2 I 5 

—2 0 ^ ^ ^ 3 


(«P-«r)v.y-/+0+«.v.y (6,63) 


The hrst term of the right-hand side can be estimated using integration by parts 


At(;0^^ 


< 


Vw ■ V (0^2) 
1/2 




< 


IV0IT 


Vw ■ Vipf^ 


+ 2 


1/2 


+ 2 


IviuHvei' 


Vw ■ Vfipf 
1/2 




1/2 


By restricting ^ G TTi, we have |V^| < |^| hence 

1/2 / /• \ 1/2 


Awtjjf^ 


< 2 


iVinpe'j yj +y (|Vn;|2 + |V0|2)e0 (6.64) 


We now estimate the remaining terms on the right-hand side of (6.63) 


7/3 7/3 

Ml —Mo 




+ 




<^7 / |M;||0|e^<C'( / w^e 


1/2 


0 "e" 


1/2 


(6.65) 


Combining the estimates (6.64)-(6.65) with (6.63) and using that infMi > cm,uj > 0 and 
M 2 > 0, we obtain 




Cm,u 


< c 


1/2 


' \Vw\^e 

(|VM;p + |V 0 p)e 


w^e 


1/2 


0 ^e" 


1/2 


( 6 . 66 ) 
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Applying Young’s inequality twice and using (6.41) of Lemma 6.4 yields 



f ’p^e+c j 

+ Vtcp -h VV’p) 



J R„M" + J (n;2 + i9')|Ver 

<-h 

f <t’V + c( 

1 Rlr^P J K + 


hence we obtain 

I {w^ + iVinl^ + ^l^^ + Rie + I {w^ + . (6.67) 

We further restrict the choice of the test function to remove the terms depending 
on w and from the right-hand side. Given C = C{M',M,uj) > 0, dehne 7 = 
min{l, {2C)~^R} > 0. First note that C Hi, so for any ^ G the estimate (6.67) 
continues to hold. In addition, |V^| < 7 |^|, hence 

j {w^- +iv«)i"++iv^7 e<c (^j ru"+ 

Rie+^j(w^+V)e) <cj Rie +1 

After re-arranging, it follows that for any ^ G 

j {w^ + I ^ \V^P\^) e<C J Rie- (6.68) 

Finally, as the equations (6.47) hold pointwise, squaring each equation and integrating 
them against yields 

j j 

j \Ai:\^e<C j {Rl + w^)e- 

Combining these estimates with (6.67), we obtain the desired result (6.44) 

j (w;2 + iVinp + \Aw\^ + <C j Rl^. □ 

Proof of Lemma 6.6. Case 1. Suppose (B) holds, so m* G M.TTk{M,uj) for some k G Nn. 
By Corollary 3.3, for i G {1,2} 

ll'*^i|lHYj(R3) R (R3) — C{k,M,L)j'). (6.69) 

Using integration by parts, we shall obtain integral estimates for derivatives of w in terms 
of derivatives of Aw. We will use the Einstein summation convention throughout this 
proof. 

To begin, we approximate w G by smooth functions Wh E C°°(M^) such that 

for all 1/91 < fc + 4, d^Wh converges to d^w pointwise [26]. This approximation is necessary 
in order to obtain estimates for d^w when la I = k + 4. 
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Fix ^ and let \(3\ = k' < k + 2. Then using integration by parts gives 



= / diid^Whdjjd^Whi 


= - j did^Whdijjd^Wh^"^ -2 j did^Whdjjd^Wh^di^ 

= / + 2 / aM,o^«.^ 3 A - 2 / 

= f + 2 f did^Whdijd^Wh^dj^ - 2 f did^Whdjjd^Wh^diC 

\ a \=2 

Summing over |/9| = k' and rearranging yields 

I Y1 ^ ^ |A9^».k|V 


\o'\=k'+2 


'\=k' 


did^Whdijd^Wh^dj^ - j did^Whdjjd^Wh^diA . 

\l3\=k'i,j=l ^' 

Then, using that ^ G C hence |V^| < |^|, we can estimate the right-hand side 
using Holder’s inequality, 

I Y, ^ f T, \^8i>w,pe 

\a\=k'+2 '' |/3|=fc' 


Idid^WhWdijd^Whl^"^ + / \did^Wh\\djjd^Wh\C 

If Y |sVkp5" + c[ / J] / J] \a^w,p( 

l^,l_t,/ 1 O 1/3. I — L,/ I/Q_ l — l,/ 1 1 


|q'|=/c '+2 


\0i\=k' 


|/32|=fc'+l 


Re-arranging this and letting h ^ 0, we obtain 


I —h./_LO ^ 1/3. I —h./ ^ 1/3^1 —II 


|q;|=/c '+2 


|/3i|=fc' 


|/32|=fc'+l 


(6.70) 


Using an identical argument, we obtain similar estimates for if, for k' < k, 

Y [\3‘-<pw<c( f Y f E 

I —ih/-L9'^ \^l\=k' \Rr.\ — hf^A ' 


(6,71) 


\a\=k'+2 


\p2\=k'+l 


In the case k' = 0, combining (6.70), (6.71) and (6.44) of Lemma 6.5 yields: there exists 
U, 7 > 0 such that for all ^ G 

f ^(|a“w ;|2 + 

| a |=2 

< U / (I Vinp + lAinp + \ViP\^ + \Aif\^) e<C [ (6-72) 
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We will now provide estimates for the right-hand terms of the form d^Aw, d^A'ip. Recall 
(6.47) 


A ^ f 7/3 7/3\ ,01+02 , U 1 +U 2 

-Aw = - ^ j + - - - W + - - - ij =: / i , 

-Atjj = An {ul - Ml) Rm =■■ /2- 

From (6.69) it follows that /i G f 2 G Let |ai| = ji < fc + 2, |« 2 | = 

J 2 + k, then differentiating (6.47) yields 

\d°‘^Aw\<C{ji,M,u) ^ [\d^^w\ + , (6.73) 

\A\<ji 

\d^^A'ilj\<C{j2,M,u) . (6.74) 

l/32|<j2 

Sqnaring (6.73)-(6.74), snmming over partial derivatives and integrating against we 
dednce 




|Ql|=jl 


\A\<3i 




|« 2 |=i 2 


\02\<j2 


Snbstitnting (6.75) into (6.70) gives for i < A; -|- 4 


/ ^ < C- / r Y Y 

\a\=ii ^|/3i|=n-l \02\=ii-2 ’ 

<c [ f Y^ Y^ 

'^|/3i|=n-l |/3i|<n-2 ^ 


|/3i|=n—1 l/3i|<n — 

Similarly, snbstitnting (6.76) into (6.71) gives for *2 < A; -|- 2 




A ia“V>l¥ < c A ( 5 ^ Y. 1 ^^ 

\a\=i2 |/3l|=l2 —1 |/32|=*2—2 

<c[( Y E {\d^^Rm\^ + \d^M"))e- 

^|/3i|=*2-1 \P2\<i2-2 ' 

Using (6.77) and (6.78), arguing by induction over Ai,A 2 simultaneously gives 

J Y (|0“w;p + |a"0p) <C J Y \d^Rm\^^^- 


| o |< fc -|-2 


(6.75) 

(6.76) 


(6.77) 


(6.78) 


To show the remaining estimate for the derivatives of w, applying (6.77) with Ai = A; -|- 
3,k + A yields the estimate (6.45) 

f( Y is”‘«’t+ E i9“"v>r)?"<cA 

^|ai|<fc+4 \a2\<k+2 ' \0\<k 
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Now fix 1 / G and choose ^{x) = e '>'1^ We will now show the lower pointwise lower 
bound for w and -0 


5^ |9“'«.(J/)P+ 5^ 

|«l|<fc +2 |Q: 2 |<fc 



where the constant C is independent of y and 7 . 

By Corollary 3.3, w G -0 G hence by the Sobolev embedding 

theorem [20] w G 0 G and 

ll^llc'=+2(Ri(i/)) < C\\w\\Hk+i(^Bi{y))^ 

||0llc7Hi(j/)) < C'||0l|H'=+2(Si(y))- 
We use these estimates to show (6.79) 


^ \ d ^^ w { y )\^ + ^ | a “ 20(|/)|2 

\ai\<k+2 \a2\<k 


^ + IIV'|lc'=.l/ 2 (spy)) 


= c 

< c 




Combining (6.45) and (6.79), we obtain the desired estimate (6.46) 

\d^^w{y)\^+ < C* f ^ |a^i?^(a:)pe- 2 ^l^-^l dx. 

|qi|<A :+2 \a2\<k |/3|<A: 

Case 2. Suppose (A) holds, then as mi G M.l 2 [M,uj), by Proposition 3.1 and (6.40), 

ll“l|lH4^if(R3) + ll'^lllHVf(K3) ^ C{M), 

+ II‘^2||j/2^.j(R3) < C{M'). 

The argument used to show (6.70) holds for k' < 2, so for ^ E Hi 



Then, as (6.75) holds with ji < 2, applying this and (6.70) for k' = 0 yields 

E f\a^'MV<c [ E + ia'’Tt) 

<c(^J J Y1 


|ai|<4' 


l/3l|<l 


\M<2 


(6.80) 
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Similarly, the argument used to show (6.71) holds for k' = 0, to give 

1 ^- 1^0 ^ ^ ^ / 


(6.81) 


|a 2 |< 2 ' 


m<i 


Finally, combining (6.80)-(6.81) and applying (6.44) from Lemma 6.5, we obtain the 
desired estimate (6.45) with fc = 0 


<c( f (|Am|" + |A^O^"+ f Y \{3^'w\" + S^'i’ne) <C [ Rie- 

\ -I I1^1 / J 


|ai|<4 


l/3i|<l 


The argument used in Case 1 holds for fc = 0 to show the desired estimate (6.46) 

\d°‘^w{y)\‘^ + \'ip{y)\‘^ - ^ f dx. 

|ai|<2 

We have now established all technical prerequisites to prove Theorems 3.4 and 3.5. 


□ 


Proof of Theorem 3.5. Applying Lemmas 6.4 - 6.6 with the assumption (B) yields the 
desired estimates (3.10)-(3.11). □ 


Proof of Theorem S.f. Case 1. Suppose spt(m 2 ) is bounded and m 2 ^ 0. We show 
assumption (A) is satished, so by applying Lemmas 6.4 - 6.6 we obtain the desired 
estimates (3.8)-(3.9). 

Since m 2 G L ,^„if (R^), it follows that m 2 G L^(R^) and since m 2 > 0 and m 2 ^ 0, it 
follows that J m 2 > 0. Then, dehne the minimisation problem 

V G > 0, / v^= m 2 > 0 

Jr3 Jr3 

which yields a unique solution (^ 2 , 02 ) to (2.4), satisfying U 2 > 0, using [31, Theorem 
7.19]. Applying Proposition 3.1, we obtain the uniform estimates 

ll“2||H;f^.j(R3) + ll02||r/2^.j(R3) < C{M'), 

Case 2. Suppose m 2 = M 2 = 02 = 0, then by dehnition (m 2 , 02 ) solve (2.6) and (A) is 
satished, so Lemmas 6.4 - 6.6 imply (3.8)~(3.9). 

Case 3. Suppose spt(m 2 ) is unbounded. By Proposition 6.3, there exists (m 2,02) 
solving (2.6) corresponding to m 2 and satisfying U 2 > 0. As we can not guarantee that 
U 2 > 0, we can not apply Lemmas 6.4 - 6.6 directly to compare (mi,0i) with (m 2,02)- 
Instead we follow the proof of Proposition 6.3 and use a thermodynamic limit argument 
to construct a sequence of functions (M 2 ,i?„, 02 ,m„) that satisfy (A) for sufficiently large 
which converges to (m 2 , 02 )- 

Let Rn t 00 and dehne m 2 ,R„ := m 2 ■ then as m 2 G L2^if(R3), m 2 > 0 and 

m 2 ^ 0, it follows that m 2 ,/j„ G L^(R^) and for sufficiently large / m 2 ,_R„ > 0. By 
Proposition 6.3, the minimisation problem 


/™(m 2 ) = inf I E™(M,m 2 ) 



rTFW 


{fn2. 


Rn 


= inf E 


pTFW, 


X, m2m 


V G > 0, 



5 


dehnes a unique solution (m 2 ,k„, 02,k„) to (2.4), satisfying M 2 ,r„ > 0 and 

||w2,R„||Ff4^.j(R3) + ||02,R„||Ff2^.j(R3) < C{M'), 


(6.82) 
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where the constant is independent of i?„. Passing to the limit in (6.82), there exist 
U 2 e 02 £ such that, respectively, along a subsequence 

converges to M 2 , 02 , weakly in and ih^(i?/j(0)), strongly in and 

L^(i?/j(0)) for all i? > 0 and for all |q!| < 2, d°‘u 2 ,R„, 02 ,converges to d°'u 2 , 02 pointwise. 
It follows that {u 2 , 02) is a solution of (2.6) corresponding to m 2 , satisfying M 2 > 0 and 
(3.7) 

l|U2||ir4^.j(R3) + ||02||//2^,^(]R3) < C{M'). 

In addition, (m'^,0'0 = (mi,0i) and (^ 2 , 02 ) = (u 2 ,m„, 02 ,m„) satisfy assumption (A) for 
large Rn, so by Lemmas 6.4 - 6.6 that there exist C, 7 > 0, independent of such that 
for large and any ^ G 

|a“i(Mi-M2 ,rJP + Y l'9“H0i-02 ,rJPV^ < C' f (ml-m2,RJ2^^ 

\ \^ A I l--^0 ' w 

|q:i|<4 |q:2|<2 

(6.83) 


and for any 1 / G M^, 



(6.84) 


Using the pointwise convergence of (m 2 ,r„, 02,r„) to (m 2 , 02), applying the Dominated Con¬ 
vergence Theorem and sending —)■ cxd in (6.83)-(6.84) we obtain the desired estimates 

(3.8)-(3.9). □ 

6.3. Proofs of Applications. The proof of Proposition 4.1 is an application of Theo¬ 
rem 3.4. 

Proof of Proposition 4-1 ■ Observe that (mi,0i) = (m, 0) and (m 2 , 02 ) = {un,(t>n) satisfy 
the conditions of Theorem 3.4, there exist C, 7 > 0, independent of D, such that for all 
y eMf 

Y \9°{u - un){y)\‘^ + 1(0 - 0 o)(i/)P f |(m - mn)(x)pe“^^l^“^l da;. 

|a|<2 

Now let 1 / G D, d = dist(i/, 90) and observe that m—m^ G L^j^;f(M^). Since sup^.^^ e- 27 hl < 
C inUgA for any A C Bi(z),z G M^, with C = C(j) independent of z, we have the 

bound 

f |(m - mf^)(a;)|2e-2^l^-^l da; < C (ll’u|li2^.j(R3) + ll”^H|li2^.^(R3)) [ da;. 

JBd{yY ^ ™ ™ ^ JBdivY 

Therefore, we obtain the desired estimate (4.1) 

f |(m — mQ)(a;)pe“^'^l'^“^l da; = f |(m — mn)(a;)pe“^'’'l®“^l da; 

Jrs Jn<^ 

< f m(x)2e-2^l^-^l da; < CM^ [ da; < CM^ [ da; 

= C]VP{1 + 

for any given 0 < 7 < 7 , where C = Cify^'y)- D 

Next, we now prove Corollary 4.2 as a direct consequence of Theorems 3.4 and 3.5. 
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Proof of Corollary 4-3. Let k E Nq and mi, m 2 G and recall the estimate 

(3.11) of Theorem 3.5, that there exists C, 7 > 0 such that 

\d°‘^w{y)\‘^ + ^ \d°'^i/j{y)\'^ < C j \d^ Rm{x)\^dx . (6.85) 

|ai|<fc+2 \a2\<k 


( 1 ) Rm having compact support is a special case of exponential decay, hence we consider 
only the case Yl\y\<k Rm{x)\^ < h-^l. it is straightforward to see that there exists 

C, 7 > 0 , depending on 7 , 7 ', such that 



■27'h-2|g 


■2'i\x-y\ 


dx < 


( 6 . 86 ) 


Hence (4.3) follows immediately from combining (6.85) and ( 6 . 86 ). 

( 2 ) Suppose that Rm satishes the algebraic decay X]|/ 3 |<fc \d^Rmix)]"^ < C{1 + |a;|)“^'’. 
It is again elementary to show that there exists C = C'(r) > 0 such that for all y 



(1 + |a;|)-2'-e-2T'l^-^l dx < ^(1 + \y\)-^^. 


(6.87) 


Combining (6.87) with (6.85) gives the desired estimate (4.4). 

(3) Now suppose that Rm E and recall (3.10) of Theorem 3.5, that there exists 

C, 7 > 0 such that for all f E 

[ ( Y1 <C [ ^\d^Rm\^^^- ( 6 . 88 ) 

^ |ai|<fc+4 \a2\<k+2 ' \l3\<k 

For any 0 < 7 < 7 , the function ^y{x) = E Hy. Then substituting ^y into ( 6 . 88 ) 

yields 

/ { \d°'^w{x)\^ + ^ dx 

'^|«l|<fc+4 \a2\<k+2 ' 

\d^Rm{x)\^dx < C f ^ Id^Rmix)]"^ da;. 

\0\<k 

Sending 7 —)■ 0 and applying the Dominated Convergence Theorem yields the desired 
estimate (4.5). 

Under the assumptions of Theorem 3.4 with k = 0, other than applying Theorem 3.4 
instead of Theorem 3.5, the proof is identical. □ 



We turn to the proofs of the charge-neutrality estimates. 

Proof of Theorem 4-3. Recall that pi 2 = mi — u\ — m 2 + u^- Let i? > 0 and choose 
Pr E C)?°(R^) satisfying 0 < < 1, = 1 on i?_R(0), (pR = d outside i?i^+i(0) and 

||<yC/j||i4/2,oo(R3) < Cy,. Let Ar := Br^i{0) \ Br{0). Recall (6.47b), that the difference 
ijj := (j)I — (j )2 solves 

— Afj = 47rpi2 (6.89) 

pointwise. Testing (6.89) with ipR and using integration by parts yields 

I P 12 PR = I 

JBr+PO) 471 Jar 

Since ipR = 1 on i?_R(0), we deduce 

I P12 = / 'ff^PR - I P12PR, 

JBr{0) 471 JAr 
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and hence 


'Hfl(O) 


Pl2 


<c (|mi - m2| + \ui - M 2 I + 101 - 02|), 
J Ar 


(6.90) 


where C depends only on c^. Observe that < CR^. 

(1) By (4.3) of Corollary 4.2 there exists C, 7 > 0 such that 

1(01 - 02)(x)| + |(mi - m2){x)\ + |(mi - U2){x)\ < 

Then using (6.90) we deduce 


'Hfl(O) 


P 12 


<C (|mi - 7712 ! + |mi - M 2 I + 101 - 020 
J Ar 

<C [ dx < 0(1 + 


(6.91) 


which implies (4.9) for any 0 < 7 < 7 . 

( 2 ) Suppose now that |(mi — m 2 ){x)\ < 0(1 + |x|)“'’, then using (6.90) we obtain 


'Br{0) 


P 12 


< C / (1 + 1x1)-" < 0(1+ i?) 


2-r 


'A. 


(3) Suppose mi — m 2 G L^(M^), then by Corollary 4.2, ui — M 2 , 0 i — 02 G 
hence by Proposition 3.1 uf — u 2 G L^(M^). Taking the Fourier transform, f{k) = 




.—2'Kik-x 


dx, of (6.89) and rearranging, we obtain 


Pi2{k) 

w 


= 7r0(fc) G L\R^). 


Arguing as in [10] we show that 0 is a Lebesgue point for pi 2 . For e > 0, 


la 


\Pi2{k)\ dk < 


f \ky dfc) ' 

Be( 0 ) / 


^BepS) |-^£(0)l 

< Oe^'^^1101 — 02 ||l2(R3), 

which tends to 0 as e —)■ 0 , as claimed. 


|Pl2(/s)P 


ae(0) 


\k\^- 


dk 


1/2 


□ 


6.4. Proof of energy locality. To prove Theorem 4.4, we hrst establish the existence, 
uniqueness and regularity of the solutions to the linearised TFW equations. 

Fix Y = {Yj)j^jq G yi, 2 [M,u) and let m = my G Ati 2 (M,ca). Let V G M^\{0}, /c G N 
and for h G [0,1] dehne 

Y’^ = {Yj + 6jkhV\j eN}, (6.92) 

and the associated nuclear conhguration 

mh{x) = m(x) + 7 (x — — hV) — rj^x — Yk). (6.93) 

Lemma 6.7. There exist M'^ojq^oj'i > 0, such that for u’ = (ci;o,Ci;(), mn G A 4 l 2 (M', cn') 
for all h G [0,1]. In particular, Y^ G 3 ^x 2 (M', a;') for all h G [0,1]. 

Proof of Lemma 6.1. Recall that mh, 7 > 0, 7 G C°°(M^) and 7 = 1, then 
sup \\mh\\L^(Bpx)) < sup \\m\\L2(Bi{x)) + ( / p{z-Yk- hvf dz 

xgR3 xgR 3 y \JBi(x) 

< M + ||7||i2(K3) =: M'. 
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inf 

xGR3 


/ rnh{z) dz > inf 

IBr{x) . 


Br{x) 


m{z) dz — 


'Br{x) 


ri{z — Yk) dz > cuqR^ — cji — 1 , 


hence for uj' = (a;o,Ci;i + 1), rrih G ,u') for all h G [0,1], as claimed. 


□ 


As nih G ,u') for all h G [0,1], by Theorem 2.1 there exists a corresponding 

gronnd state {uh,(t>h)- Also, let (n, 0) = (mo)0o)- We now use Corollary 4.2 to compare 
{uh, 4>h) with (u, 0) and rigorously linearise the TFW equations. 

Lemma 6 . 8 . Let Y G yL^{M,u) and let m = my G JYli 2 {M,uj). Also, let fc G N, 
C G \ {0} and ho = min{l, |ld|“^}. For h G [0, ho] define 

mh{x) = m{x) + T]{x - Yk - hV) — r]{x - Yk). 

There exist C = C{M',uj'), 70 = 7 o(M',a;') > 0, independent of h and \V\, such that 

- m)(2 :)| + |(9“(0h - 0)(2:)|) + \{mh - m){x)\ < (6.94) 

l « l <2 

\\uh — m||h4(k3) + 110/1 — 0||//2(r3) < C\\mh — m||2,2(]K3) < Ch. (6.95) 

Moreover, the limits 

- r '^h-U - 4>h-4> — r 

h-^-o h h^o h h^o h 

exist and are the unique solution to the linearised TFW equations 

— Au + — 0 ^ M — M 0 = 0 , 

— A0 = 471 (m — 2uu). 

Moreover, u G i/^(M^),0 G if^(M^),m G C^(M^) and satisfy 

(|5"m(2:)| + \d’^fi{x)\) + \m{x)\ < 

|«|<2 

||'?i||h'»(r 3) + ||0||iy2(K3) < C'||m||i2(R3). 


(6.96a) 

(6.96b) 

(6.97) 

(6.98) 


Proof of Lemma 6.8. By Proposition 3.1 and Proposition 3.2, for h G [0, ho] the ground 
state {uhyfih) satishes 


.7R3) + I|0/i||h2 .7R3) < C{M'), 


inf Uhix) > cm',u' > 0 , 

a;gR3 

independently of h. From (6.93), it follows that 
\{mh - m){x)\ = \7]{x -Yk - hV) - r]{x - Yk)\ 


(6.99) 

( 6 . 100 ) 


<h\V\ \V7]{x-Yk-thV)\dt< \V7]{x-Yk-thV)\dt. (6.101) 
Jo Jo 

For all h G [0,ho], spt(m/i — m) C Bjig+fiYk), so by Corollary 4.2 and (6.101) it follows 
that there exists 70 > 0 such that 

Y l'9“(Mh -u){x)\ + 1 ( 0 /, - 0 )(a;)| + |(m/, - m){x)\ < 
l«l <2 


( 6 . 102 ) 
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and (6.95) holds 

\\uh - m||r4(ir3) + \\(j)h - 0||h2(r3) < C\\mh - m||L2(iR3) < Ch. (6.103) 

Due to the uniform estimates (6.99)-(6.100) and (6.101), the constants appearing on the 
right-hand side are independent of h. 

We now show 

^ \dP{(j)h - (t)){x)\ < (6.104) 

l«l<2 

Observe that for h G (0, ho] as spt(m/i — m) C by the triangle inequality 

X e implies B 2 {x) C 5|j^+i(Yfe). Consequently, for x G 5|j^+3(Dfc) 

\\mh - "i|lcFi/2(B2(x)) = 0, (6.105) 

and for x G Bn^^siYk), by (6.93) it follows that 

W'^h — ’^||c0’1/2(£2(x)) ^ 2||?7||c'0,l/2(£2(a:))- (6.106) 

By (6.105)-(6.106) we deduce that x i—)■ \\mh — mo\\(jo,i/ 2 (^B 2 {x)) is a bounded function with 
support in iJ^g+oCW), hence there exists C > 0 such that 

W^h - m\\co,i/2^B2ix)) < (6.107) 

Then we apply the Schauder estimates [26, Theorem 10.2.1, Lemma 10.1.1] together with 
(6.102) and (6.107) to estimate 

110/x - 0|lc'2a/2(Bp3,)) < C {\\mh -m-u\ + U^\\cO:i/^(B 2 ix)) + \\4>h - 4>\\L'2iB2ix))) , 

< c {\\mh - m\\c0.i/2(^B2ix)) + Il'Wft ~ u^\\co^/^iB2ix)) + W^Ph - 4 ’\\l^(B 2 {x))) , 

< C [\\{Uh + u){Uh — u)\\cO,i/2(^B2(x)) + ^ ) 

< C {\\Uh + u\\cO,i/2^B2{x))\\'^h — u\\c°’^/'^{B2{x)) ^ ^*’0 • (6-108) 

Applying the Sobolev embedding C^'^^‘^{B 2 {x)) ^ H‘^{B 2 {x)) and using (6.100), it follows 
that 

Il'Wft + '^IIc'0.i/2(B2(x)) ^ C^Uh + u\^m{B2{x)) < C (^||'*I/i||h2^.j(R3) + ll'^llH, L if(«^)) — 

(6.109) 

Applying (6.109) and (6.102) to (6.108), we obtain the desired estimate (6.104): for any 
multi-index a satisfying |a| < 2 

|^“(0/1 - (p){x)\ < 110/1 - 0|1iv2.°°(Ri(x)) < 110/1 - <P\\c^^/2(Bi{x)) 

< c {\\Uh + M|lc 0 .i/ 2 (B 2 (a,))||n/i - m||co.i/ 2 (B 2 (x)) + 

< C {\\uh - n||wLoc(S 2 (.)) + (6.110) 

We will show next that there exist u G i7^(M^),0 G i7^(M^) such that 
converge to u, 0 respectively, weakly in and strongly in H^{Bb{0)) and 

H^{Bb{0)) for all i? > 0 and pointwise almost everywhere, along with their derivatives 
as h ^ 0. 

First consider any decreasing sequence h„ —)■ 0, then there exists a subsequence (still 
denoted by h„) such that converge to u G i7^(M^),0 G i7^(M^) respectively, 

weakly in and strongly in H^{Bb{0)) and H^{Bb,{0)) for all 72 > 0 and 

pointwise almost everywhere, along with their derivatives. In addition, it follows that 
(u, 0) satisfy (6.97)-(6.98). 
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We now verify that the limiting fnnctions are independent of the choice of seqnence. 
First, observe that by passing to the limit as —)■ 0 in the eqnations 


» I Uh„ — U \ “h 

_ A ( — - H-— 

hr,. J 3 h, 




= 0 , 




hr 


hr 


it follows that (m, 0) solve the linearised TFW eqnations (6.96) pointwise, 

'35 


-Am 


9 


-u 


4/3 


M — M(/) = 0 , 


— A 0 = dvr (m — 2 mm) , 

where m[x) = lim -;-= ~\/r][x — Yk) • V. 

hn —^0 hyi 

Clearly m is independent of the seqnence hn- Applying [7, Corollary 2.3], it follows that 
the (m, 0) is the nnique solntion to the linearised system (6.96), hence is independent of 
the seqnence (h„). It then follows that converge to m, 0 as h —)■ 0 as stated 

above. □ 

We are now in a position to prove Theorem 4.4. 

Proof of Theorem 4-4 ■ We will repeatedly nse the fact that there exists C, 7 > 0 snch 
that, for all h G [ 0 , ho], p G [ 1 , 2 ], 


(1 + m;,(x) + |V 0 ,,(a:)|)^’e-^“l"-^'=le-^l"-^^l da; < 


( 6 . 111 ) 


which is a conseqnence of the nniform bonnds on rrih, (ph and of ( 6 . 86 ). 

Fnrther, we reqnire that there exist C, 7 > 0 snch that, for j G M, h G (0, ho], x G M^, 

ipj{Y^]x) - (Pj(Y;x) 


h 




( 6 . 112 ) 


which follows directly from (4.19c). 

For z = 1, 2 and j G Id, consider the difference 

E}(Y^)-E](Y) f g,{Y^-x)ip,{Y^-x) - S,{Y-x)if,{Y-x) 


h 


h 


dx 


S,{Y>^-,x)-Si{Y-,x) 

h 


(Pj{Y^] x) dx 


+ ('MAApAIih]) dx 

We wish to show that the limit of (6.113) exists as h —?■ 0 to obtain 


dE) 

w: 


^{Y;x)^j{Y;x) dx + [ gi{Y;x)^{Y;x) dx, 
Oik Jr 3 oik 


(6.113) 


(6.114) 


^^(F; ■) = 2Vm • Vm + ^M^/^M + - u^) + ^0(m - 2 mm), 

uYk o 1 1 

|^(F; ■) = 2Vm ■ Vm + ^m'/^m + ^V0 • V0. 
oYk 3 dvr 


(6.115) 

(6.116) 


where 
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(6.117) 


h 


= V{uh + m) ■ V 


Uh-U 


,10/3 _ 10/3' 


Uh' -u 


+ TT - ul) - -0(m - M^) 


2h 

= V {uh + u) ■ V 
, ^ f (ph-(p 


Uh-U 

h 


,10/3 _ 10/3' 


% 


^ Am-u^) + -(t)h 


h 

rrih — m — u\ + 

h 


(6.118) 


It follows from (6.118) and pointwise convergence of Uh, Y^Uh, (ph to u, Vm, 0 and 


Uh-U 

h 


’ V (^) , to M, VF, 0, m as h ^ 0, that (6.115) holds 


, 8i(Y'^;-) -8i(Y;-) _ 10 1-, 2^ 1,/— -x d8i 

hm- - -= 2Vm ■ Vm + —vJ'^u + -0(m — u^) + -0(m — 2uu) = ——. 

>^0 o 2i 2i k 

Applying (6.94) to (6.118) yields 

\8i{Y^]x) - 8i{Y]x) \ < C {\{uh - u){x) \ + \ V{uh - u){x) \ + \ {mh - m){x)\) 

+ C(l + m(a;))|(0/,-0)(x)| 

< Ch{l + m(a;))e“'^°l^"^'=L (6.119) 

Combining (6.119) and (4.19b), we deduce 


-fiiiy-, 


Xl 


< C(1 + m(x))e 


-7oh-Yfc| „-7h-Yj| 


( 6 . 120 ) 


hence by (6.111) and the Dominated Convergence Theorem, 


JR3 OYk 

It follows from (6.120) and (6.111) that 

f d8i 


8,iY'^;x)-8iiY;x) 

h 


^ (Pj(Y;x) dx. (6.121) 


dY, 


{Y ; x)ipj{Y ; x) dx 


<C (1 + m(x))e-^“l^-^'=le-'^l^-^^l dx < 

Jms 

( 6 . 122 ) 


It remains to show that (6.113) converges using (6.111) and (6.112). As (pj(Y;x) is 
differentiable with respect to 10, for all x G 


£i(y;i)^(y;a:) = limfi(y 

uY k /i—lO 




X 




X 




and combining (6.117) with (6.112) implies 

£.(y;,) I < C(1 + 


-7o|x-Y'fe|g-7|x-yj| 
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hence by (6.111) and the Dominated Convergence Theorem, 


Ei{Y]x)^^{Y]x) dx = lim / 8i{Y]x) 
oYk 


dx, 


and 


Si(Y-,x)^{Y-,x) ,lx 
OYk 




Combining (6.122) and (6.123) yields the desired estimate (4.21). 
The second case, using S 2 instead of Si, is analogous. 


Proof of (4.22). We will use 


JGN 


,-7lD-n| 


< 00 , 


(6.123) 


□ 


(6.124) 


which is a consequence of (HI) and that Y G yL 2 [M,uj). Then for i G {1,2} 

dS, 


Y 

jGN 


dE) 


OYk 


JGN 

jeN 


{Y; x)ifj{Y ; x) dx 


dYk 

e < 00 , 


isN 


e,{Y;x)^^{Y-,x) di 


hence by the Monotone Convergence Theorem, the sum is well-dehned 


vEh 

^ dYk 




x) ) dx. 


As satishes (4.19a) for all h G [0, ho], it follows that 


JGN 


and consequently. 


^ dYk 
jen 


dS, 

dYk 


(Y ; x) dx. 


Now consider the difference of (6.115)-(6.116) 




(6.125) 


and applying integration by parts yields 

' dSi dS2 


dYk dYk 


(Y ; x) dx = 


-<p{m — u^) + -(t){m — 2uu) — • V0 


1 

Stt 

1 

Sir 


A0) + 0(-A0) - 2V0 • V0) 
(2V0 • — 2V(j) ■ V0) = 0. 


In addition, since 


— / V0 • = ^ / 0(—A0) = / 0(m — 2uu) 


Air 


Air 
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and since u solves (2.6a), —Au + — 0m = 0, the desired resnlt (4.22) holds: 

f ^ [ ( Vm • Vm + — 0MM ] + [ (j)m = f (j)rn. □ 

Jr 3 cYfc Jr\ 3 / Jr3 Jr3 

Finally, we establish (4.16)-(4.17): if the partition functions (fj are invariant under 
permutations and isometries, then so are the site energies. 


Lemma 6.9. If the partition is permutation and isometry invariant (4.23)-(4.24), 

then for i = 1,2, for any bijection P : N —?■ N, isometry A : —)■ j E N and 

YeyL^{M,cj) 

E}{YoP) = E}{Y), (6.126) 

E}{AY) = E}{Y). (6.127) 

Proof of Lemma 6.9. Let Y G yL 2 (M,u) and m = my, then as P : N —)■ N is a bijection, 

mY^p{x) = '^r]{x - Yp.) = '^r]{x - Yj) = my(a;). 

jeN jgn 

Since (2.6) has a unique solution, {uy^cfy) = (myop,0yop)- Consequently, the energy 
densities agree, £i{Y o P; ■) = £i{Y] •). Together with (4.23) this implies (6.126). 

We now show isometry invariance (6.127). First consider a translation Ai{x) = x + c, 
for c G M^, then 

LnAiy(x) = ''^^r]{x — Yj — c) = mY{x — c) = mY{Af^{x)). 

jeN 

Then, by the uniqueness of the TFW equations, it follows that {ua^y, 0Aiy)(-) = {uy, 4>y){,- — 
c), so £i{AiY ; •) = £i{Y ; • — c) and thus 


pj(7iiy) 


£i(AiY; x)ipj(AiY; x) dx 


Si(Y ; X — c)(pj(Y ]x — c) dx 


£,{Y-z)^,{Y-z) dz = E]{Y). 


Similarly, for a rotation A 2 {x) = Pt, R G 0(3), since we assumed that rj is radially 
symmetric. 


mA2Y{x) = ''^ri{x — RYj) = ''^r]{R{R^x — Yj)) = '^/q{R^x — Yj) = my^R^x). 

jgN jgN jgN 

(6.128) 


As (wy, 0y) solve (2.6) 

— Auy H Uy — (fyUy = 0 , 

3 

— A0y = 47r(my — My), 

then by (6.128) and as the Laplacian is invariant under rotations, it follows that 
(m, 0) = (My, 0y) O Af^ SolveS 

5 

—Am + -m’^'^^ — 0m = 0, 
o 

— A0 = 47r(my o R^ — u^) = 4:7r{mA2Y — m^), 
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hence the uniqueness of (2.6) implies {ua 2 Yi (j>A 2 Y) = {uy,^y) ° ^ 2 ^- If follows that 
Ei[A 2 Y ; •) = EiiY ; hence as det(i?) = 1, a change of variables shows 


Ei{A 2 Y)= / 8 i{A 2 Y-x)ipj{A 2 Y-,x) dx= / 8 i{Y-R^x)^,{Y-R^x)dx 
JR3 JR3 

= [ 8 ,{Y-,z)ipj{Y-,z)\det{R)\dz= f 8 i{Y-, z)ip,{Y-, z) dz = EUY). 


As the site energies are invariant under both translations and rotations, they are invariant 
under all isometries of □ 
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